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Abstract. We describe here the new concept of e-Homomorphisms of Prob- 
abihstic Regulatory Gene Networks (PRN). The e-homomorphisms are special 
mappings between two probabilistic networks, that consider the algebraic ac- 
tion of the iteration of functions and the probabilistic dynamic of the two 
networks. It is proved here that the class of PRN, together with the homo- 
morphisms, form a category with products and coproducts. Projections are 
special homomorphisms, induced by invariant subnetworks. Here, it is proved 
that an e-homomorphism for < e < 1 produces simultaneous Markov Chains 
in both networks, that permit to introduce the concepts of e-isomorphism of 
Markov Chains, and similar networks. 

Introduction 

We can understand the complex interactions of genes using simplified models, 
such as discrete or continuous models of genes. Developing computational tools 
permits description of gene functions and understanding the mechanism of regula- 
tion m ■ We focus our attention in the discrete structure of genetic regulatory 
networks instead of continuous models. Probabilistic Gene Regulatory Network 
(PRN) is a natural generalization of the Probabilistic Boolean Network (PEN) 
model introduced in 0, and Pp. This model have n functions defined over a finite 
set X to itself, with probabilities assigned to these functions. We present here the 
ideas of e-similar networks, and isomorphism of Markov Chain, e-homomorphisms 
are used to describe subnetworks and similar networks, because they transform the 
discrete structure of one network to another, and the probability distributions of 
the networks are enough close, using a preestablished < e < 1 as a distance 
between the probabilities. 

1. Preliminaries 

Probabilistic Regulatory Networks A Probabilistic Gene Regulatory Network 
(PRN) (or a Probabilistic Dynamical Systems) p] is a triple X = {X, F, C) where X 
is a finite set and F = {/i, ...,/„} is a set of functions from X into itself, with a list 
C — (ci, . . . , c„) of selection probabilities, where Ci = p{fi), P We associate with 
each PRN a weighted digraph, whose vertices are the elements of X, and ii u,v ^ X , 
there is an arrow going from u to v for each function fi such that fi{u) = v, and the 
probability c; is assigned to this arrow. This weighted digraph will be called the 
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state space of X. In this paper, we use the notation PRN for one or more networks. 
li X = Xi X ■ ■ ■ Xn is the product of n sets of possible values of the variables, 
then with the vector function / = (/^, • ■ • , /") we associate a digraph F, called 
dependency graph, with vertex set {1, . . . ,n}. There is a directed edge from i to 
i if xi appears in the component function /-'. For a PRN, we have a dependency 
graph (dep-graph) for each function, then we superpose all the dep-graph and that 
is the low level digraph of our PRN jfj 



Example. Suppose we have two genes with two values that we denote as usual 
{0, 1}, that is this PRN is a very simple PEN. The set of boolean functions F is 
the following: 

F = {h{x\.x2) = (xl,Q)j2{xl,x2) = (l,x2), 
f3ixl,x2) = (l,0),/4(xl,x2) = {xlx2,x2)}, 

and the probabilities are {.21, .22, .34, .23}. Therefore, the PEN X = {X, F, C) has 
the following state space, dependency graph, and transition matrix. 



.66 

O (0,0)^ 

.34 ^^.23 

O (1,0)^ (1,1)0-45 
State space 
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e-Homomorphisms of PRN. If C is a set of selection probabilities we denote by x 
the characteristic function over C. That is x ■ CU{0} {0, 1} such that x(c) = 1, if 
c ^ and x(0) = 0. Let - (Xi, F - (/OLi, C) and X2 = (X2, G = {gj^JLi^D) 
be two PRN. A map : Xi X2 is an e-homomorphism from Xi to A2, if for 
a fixed real number < e < 1, and for all fi there exists a g^, such that for all u, v 
in Xl, 

(1) (/) o = gij o (/); (2) maxu,v\cfi{u,v) - dg^{(j>{u), (f)iv))\ < e, and 
(3) xK, (</-(«),</-(«))) >x(c/,Ki^)). 
If (j> : Xl ^ X2 is a bijective map, and dg-{4>{u),4>(v)) — Cf.{u,v), for all /i, gj, u, 
and V in Xi\ then cj) is an isomorphism. 

If we denote by p{u, v) =Y,f, Cf.{u, v) and p{(j){u) , (j>{v)) = J^g^ dg^ ((^(m), (/)(«)), 
then condition (2) implies that \p{u,v) — p(0(w), </)(?;) )| < fee, where k is the max- 
imum number of functions going from one state to another in the network. So, if 
Ti denote the transition matrix of Xi, and the entry (u, v) of Ti is p{u, v) then the 
third condition implies that: maXu,v\{Ti)u^v — (T2)^(u),4>{v) \ ^ for all possible u 
and V in Xi. 



2. Isomorphism of Markov Chains, e-Similar Networks 

Two PRN are e-similar if there exists a bijective homomorphism (j) between them, 
such that (t)~^ is also an homomorphism. Observe that </> and (/)~^ have the same 
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e. When two PRN are e-similar, the two transition matrices have the a similar 
distribution of probabiHties. 



Theorem 2.1. If cf) : Xi —> X2, and (f) are bijective e-homomorphisms, then 

max\{cfm [u, 
for all m > 2; u, v, in Xi . 



Proof. If x{cf{u, f{u))) = 1, then x{dg{4'{u), (j){f{u)))) — 1, because and are 
bijective homomorphisms. By definition of e-homomorphism, g{(j3{u)) = 4>{f(u)). 
Then for m = 2, and by the Chapman-Kolmogorov equation 0, we have the 
following: 

\cp{u,.f{u))-dg.{c^{u),g\<P{u)))\ = 
\cfiuJiu))cfifiu)j\u))-dg{q^iu),gi<Piu)))dMiu)),g^m^^^ 
\cf{u, f{u))cf{f{u), f\u)) - d,{<j,{u),^{f{u)))d,{<j,{f{u)),4>{f^{u)))\ < 
By condition (2) in definition of homomorphism, we have 

< \cf{f{u), f^u))\e + \dM^)^ m^Me < 2e. 

Then we proved that \cf2{u, f^{u)) — dg2{(j){u), g'^{(j){u)))\ < 2e. 
Using this property, and mathematical induction over m, we can conclude that 
our claim holds. □ 

Corollary 2.2. If (j) : Xi ^ X2, and (f>~^ are bijective e-homomorphisms, then the 
transition matrices Ti and T2 satisfy the condition: 

1. x{Tr)u,v = x{T^)u,v, 

2. Y.tiiiTr)u..~{W)u.,) = 0, 
for all m, u — and u = (piu). 

An e-homomorphism between two PRN determines a correspondence between 
the Markov Chains of these two networks. Here, we introduce the concept of two 
similar Time Discrete Markov Chain (TDMC). 

Definition 2.3. Two TDMC of the same size n x n: {Ti, Tf, Tf, . . .}, and 

{T2, T2 , T2 , . . .} are e-similar or e-isomorphic if there exists on e G R small 



I J y n X n 



satisfies that 



enough, such that T^' 

(1) \Uj\ 

(2) x{Ti^)ij = x{T2^)ij, for all m, where x the characteristic function. 
That is, these two TDMC simulated the dynamic of two e-similar networks. 



< e, and *y = 



Example 2.4. 

The networks with dynamic Ti and T2 are .005-similar. In fact 
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As a consequence, we obtain max\{Ti)ij — (T2)ij\ < .005, and both dynamics are 
.005-isomorphic. The steady state of Ti is tti = (0, .01926, 0, .98074), and the 
steady state of T2 is tt2 = (0, .01632,0, .98368), |9j. We can see that Itti - tts] = 
'maXi\Tri{i) — TT^{i)\ < .004. Additionally, we have 



-.001467 -.00136 

.00199 

-.00232 -.00019 

.002639 



.00006 .00277 

-.00199 

.00295 -.00243 

-.00263 



therefore max\{T^)ij ~ (r|)y | < .003. 



- T2' = 
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.000161 
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.002525 
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.002843 



.00011 


.00213 




.00073 
-.00253 
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.002843 



and max \(Tl)^J - (If)ijl < .004. In the above example, the TDMC generated by 
T and T2 are .005-similar, and the networks simulated by them are .005-similar. 

3. The category of Probabilistic Regulatory Networks, and 
mathematical background 

For a e G R small enough, we have the following theorem. 

Theorem 3.1. // 0i : X^X2, and (/)2 : A2 — > A3 are ei-homomorphisms, for i — 
1, 2. Then cj) — (t>2°4'i '■ '^1 ~^ is an e-homomorphism. Therefore the Probabilistic 
Regulatory Networks with the e-homomorphisms of PRN form the category PRN. 

Proof. The Probabilistic Regulatory Networks with the PRN homomorphisms is a 
category if: the composition is an homomorphism, and satisfy the associativity law; 
and there exists an identity homomorphism for each PRN. 

(1) Let (f)i : Xi X2 be an ei-homomorphism, and let (f)2 ■ X2 X^ be an 
e2-homomorphism. If /i G ^3, g S and f € Xi are functions in each PRN, and 
such that 4>io f = go 4>i and (1)2 ° g = ho (f>2, then we will prove that: (jj o f = ho (p. 
In fact, 

{(f)2 o (f>i) o f = cj)2 o [cj)i o f) ^ cj)2o{go cj)i) = [cj)2 o g) o 01 = 

{h O 02) O 01 = /l O (02 o 0l). 

(2) To verify the second condition for e-homomorphism, we do the following. If 
Cf{(j){u), 0(w)) 7^ 1, with u, V = f{u) £ Xi, for some / S <Y2,then we will prove that 
there exists an e < 1 such that 

\cf{u,v) -t,i(0(u),0(i;))| < e. 

by part (1). We denote by u = </'i(w), v — </'i(w)- 

\cf(u,v) - dg(0l(u),0l(u)) +dg{(j>l{u),(f>l(v)) - ift(02(u),02(l)))| < 
\Cf{u,v) - dg0i(u),0i(w)| + |dg(0i(u),0i(w)) - th{<j>2{u),<l>2{v))\ < 

Therefore our claim holds, \cf{u,v) — th{4'{u),(t){v))\ <ei+e2- 

(3) We want to prove that x{th{4>{u), 0(f))) > x(c/(u, v)). Suppose that x(c/(u, v)) 
1. Then, since 0i is an homomorphism of PRN, we have that 

X(dg(0l(u),0l(w))) > X(C/(U,W)) = 1 
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Since (j)2 is an homomorphism of PRN, wc obtain that 

X{th[(^{u), (/.(«))) = x{th{MMu)),c^2{'t^i{v)))) > xicfiM^), (0i(«)) - 1- 

Therefore we have that x(i/i (02 ('/'i(u)), 02 (0i («)))) = 1- 

Then the composition of two PRN-homomorphisms is an homomorphism. 

The associativity and identity laws are easily checked, then our claim holds, and 
PRN is a category. □ 

For proofs of the following theorems see [2] 

Theorem 3.2. LetXixX2 = {XixX2,H,E) be a product of PRN Xi = (Xi,F,C) 
and X2 — {X2,G, D). If Si : X ^ Xi are two PRN-homomorphisms, then there 
exists an homomorphism S : X Xi x X2, such that (pi o 5 = Si for i = 1,2. That 
is, the following diagram commutes 

Xi X X2 

01 s 4>2 

/ T \ 



X^'^X^X2 



This homomorphism is unique. 



Theorem 3.3. Let Xi®X2 = (Xi xXa, F, E) he a product of PRN Xi = (Xi,F, C) 
and X2 = {X2,G, D). If ji : Xi ^ X are two PRN-homomorphisms, then there 
exists an homomorphism 7 : Xi © X2 X , such that j o i-i — ^i for i — 1,2. That 
is, the following diagram commutes 

Xi®X2 

/ i \ 



This homomorphism is unique. 



Xi ^ X ^ X2 



4. Subnetworks 



A subnetwork Y (- X oi X ~ {X, F, C) is an invariant subnetvirork or a sub- 

PRN of X if f,{u) e Y for aU u e Y, and /, e F. Sub-PRNs are sections of a PRN, 
where there aren't arrows going out. The complete network X, and any cyclic state 
with probability 1, are sub-PRNs. An invariant subnetwork is irreducible if doesn't 
have a proper invariant subnetwork. An endomorphism is a projection if n'^ — it. 

Theorem 4.1. // there exists a projection from X to a subnetwork y then y is an 
invariant subnetwork of X . 

Proof. Suppose that there exists a projection tt : X ^ Y . If y G y, by definition 
of projection iT{y) = y, and fi(Tr{y)) = Tr{gj{y)). Therefore all arrows in the 
subnetwork Y are going inside Y, and the network is invariant. □ 

4.1. Constructing a PRN with real data. Here we developed a method to 
construct a PRN. In this case, we suppose that the information given by the ex- 
periment is a dependency graph and a time series data, see Figure 1, and Table 1. 

Additionally, we know that this information is noisy, and the first gene has three 
values, meanwhile the other two genes take only two {0,1}, so X — {0,1,2} x 
{0,1}'. 
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Table 1. Time series data 



To determine the partially defined functions: /i, /2, /a over the finite field with 
3 elements Z3, we use the algorithm introduced in That is: the first variable 
xl G Z3, meanwhile the other two genes X2, and 2:3 are in Z2. 

For example with the first function /i = (/n, /12, /13) we do the following. We 
represent the functions with polynomials over the variables given by the dependency 
graph, and the operations + and • are the usual in the finite field Z3. Then, the 
second component function 

/i2(a;l, x2, x3) = a + bxl + cx2 + dx3 + exlx2 + gxlx3 + hx2x3 + txlx2x3 



takes the following table of values. 



/12 (mod 2) 


1 
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x2 


1 











x3 
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Evaluating, we 



obtain the following linear system, where "=" means congruence (mod 2): 



a 

a - 



-2b + c + 2e 
2b + d + 2g 
a + 25 



Then reducing modulo 2, we have a — c — d = 0, and 6, e, g, h, t, are free 
variables. So, one of the solution is /12 ~ xl{l + a:2 + a;3 + x2x3), (mod 2). Using 
this method, we obtain the following functions: 

fi{xl, x2, x3) = (a;i, a;l(l + x2 + x3 + x2x3) , X2) , 
f2ixi,X2,X3) = (a;l,a;3,0), 
/3(a:i, 3^2,2:3) = (xla:2, a;2, a;3), 

and they have the probabilities ci — .23, C2 — .34, C3 = .43. 

The state space of ^? = {X,F,C) is in Figure 1. The network has 12 states. The 
only fixed point is (0, 0, 0), and the state space has two subnetworks of 8 elements 
and one subnetwork of 4 elements. For each subnetwork we must have a projection. 
That is, an e-homomorphism tt : X Y , must exist for each subnetwork Y. That 
is, the converse of the Theorem 14. II could be true in some cases or with some little 
changes. 

In particular, for the sub-PRN 3^1 = {Yi; F; C} with 

Y, = {(1,0,0), (1,0, 1), (1, 1,0), (1, 1, 1), (0,0,0), (0,0, 1), (0, 1,0), (0, 1, 1)}, 

a projection tti : X ^ Yi exists, in fact: 7ri(a;l, a;2, a;3) = {xl, x2^ x3) if xl = 0,1; 
and 7ri(a;l, a;2, a;3) = (0,a;2,a;3) if xl — 2. With this projection, it is possible to 
consider the first gene with only two values: {0, 1}. 
For the sub-PRN 3^2 = {Y2;F; C} with 



Y2 = {(2,0,0), (2,0, 1), (2, 1,0), (2, 1, 1), (0,0,0), (0,0, 1), (0, 1,0), (0, 1, 1)}, 
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Figure 1 . State Space of G 



a projection tt2 : X ^ Y2 doesn't exist, because the first function /i. So, taking a 
subnetwork of the whole PRN but without the function /i and a new assignation 
of probabihties we have a new PRN Q = {X; f2, fs; d2, d^} and a projection ^2 : 
Q ^ y2 exists, and it is given by: 7r2(a;l, x2, x3) — (xl,a;2,x3) if xl = 0,2; and 
7r2(a;l, x2, xZ) = (0, x2, a;3) if xl = 1, where 3^2 — {^2; /2, h] ^2, ^3} The projections 
■Ki are .57-honiomorphisms. These two subnetworks and 3^2 are not similar. 



4.2. Future work. The construction of a mathematical model for the genetic reg- 
ulatory network ENDOMESODERM GENE NETWORK, described in HE], will 
be developed in the future. The subnetwork, "Mat-Act" is formed with the action 
of two genes called: Mat-c/3 and Mat-Otx over eight genes: ECNS, GSK3, Wnt8, 
N/3-TCF, Bliml-Krox, Nucl, KRL, Pmarl; whose interaction is during 21 hours. We 
will use the above methodology, for the genetic network with the dependency graph 
in Figure 2, obtained in Biotapestry 
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Figure 2. Mat- Act Dependency graph 
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